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AB STR ACT 

The itraii wan (actor (SWF) signal, which is ths output of an 
ultrasonic tasting system whsra the transmitting and receiving transducers are 
coupled to the same (ace of the test structure, is analysed in the frequency 
domain. The SWF signal generated in an isotropic elastic plate is modelled as 
the superposition of successive reflections . The reflection which is generated 
by the stress waves which travel p times as a longitudinal (P) wave and s 
times as a shear <S) wave through the plate while reflecting back and forth 
between the bottom and top faces of the plate is designated as the reflection 
with p, s Short- time portions of the SWF signal are considered for obtaining 
spectral information on individual reflections. If the significant reflections 
are not overlapped, the short- time Fourier analysis is used. 

A summary of the relevant points of homomorphic signal processing, which 
is also called cepstrum analysis, is given. Homomorphic signal processing is 
applied to short- ti'je SWF signals to obtain estimates of the log spectra of 
individual reflections for cases in which the reflections are overlapped. 

Two typical SWF signals generated in aluminum plates are analysed. In 
Test I, the SWF signal contains nonoverlapping reflections. In Test II, the 
thickness of the plate is approximately half that in Test I, and the 
t elections are overlapped Multiple reflections of only F waves are 
considered It is observed that the log spectra of two corresponding 
reflections with the same ray angles and approxiamtely the same total ray path 
lengths in Tests I and II have approximately the same characteristics except 
for their overall scales, as the theory of wave propagation in isotropic 
elastic plates suggests From the differences between the ovetall scales and 
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the values of p of the two corresponding reflections in Tests I and II, 
experimental reflection coefficients for P wavs to P wavs reflections are 
obtained for the different ray angles considered. They show good agreement 
with the theoretical reflection coefficients for plane P waves. Also, it is 
observed that the frequency response at the maximum value of the log spectrum 
moves toward higher frequencies in the reflections with larger p, namely with 
smaller ray angles This is consistent with the behavior of the directivity 
function . 

The potential use of these results for the nondestructive evaluation of 
plates containing cracks perpendicular to their faces is demonstrated in a 
specific test 
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1 INTRODUCTION 

Tha measurements of nondestructive evaluation <NDE) parameters in an 
ultrasonic testing <UT> system may be conducted directly on tha input and 
output signals or indirectly on a processed signal Cl i . These indirect 
techniques usually take the form of a transformation of a signal into another 
signal on which the interpretations of tha measurements for NOE may be 
easier The Fourier transform is one of the most commonly used signal 
processing analytical tools. 

Recently, UT systems in general have tha input and output signals in the 
form of discrete- time signals (sequences) and they are processed on a digital 
compute.'. This is faciliated by the availability of high frequency analog to 
digital converters, and the efficient algorithm for computation of the Fourier 
transform known as the fast Fourier transform (FFT) C23 . There are 

increasingly sophisticated signal processing algorithms developed to exploit 
the flexibility of the digital computer C31 Ultrasonic NOE applications of 
these algorithms are current research areas An example of a class of 
algorithms of this type is tha set of lechniques referred to as cepstrum 
analysis or homomorphic signal processing C33 , which is used in this study 

In general, a pulse having a broadband freguency spectrum or a tone 
burst is used as the input signal m a UT system The output signal may 
consist of superposed pulses or ton;* bursts, respectively, which arise due to 
reflections frt,m the boundaries of the test structure. In conventional UT 
techniques C43 such as tha pulse- echo or through- thickness transmission 
methods, the measurements are conducted directly in the time domain. Peak 


values of pulses or magnitudes of tone bursts and their time delays are 


measured. These methods are limited to output signals which contain 
nonoverlapping pulses or tone bursts 

In conventional ultrasonic spectral techniques C51, measurements on the 
output signal are conducted in the frequency domain. The magnitude or phase 
spectra of the individual pulses are analysed. These methods are also limited 
to output signals with nonoverlapping pulses. 

Recently - ry et al C4.73 introduced an ultrasonic NDE parameter called 
the "stress wave factor" (SVF) Separate transmitting and receiving 
transducers are coupled to the same face of the test structure, which is 
called the SVF configuration. An input pulse having a broadband frequency 
spectrum is applied to the transmitting transducer and the number of 
oscillations in the output signal at the receiving transducer esceeding a 
preselected voltage threshold (or some modifications thereof C8]> is defined 
as the SVF The SVF is also valid for the analysis of the output signals with 
overlapping nchoes 

The SVF analysis has shown encouraging results in several NDE 

applications C6— 93 , and its use in new applications is currently under study 
C93 This work is part of an overall effort to develop quantitative analyses 
of the SVF and computer- aided nondestructive evaluation (CANDE) capabilities. 
The stress wave transmission chracteristics of an isotropic elastic plate 
through which the transmitting and receiving transducers communicate in the 
SVF configuration are studied in CIO] and Cl 13 . The SVF signal is modelled as 
the superposition of the signals due to the stress waves which multiply 
reflect between the top and bottom faces of the plate. The esperimental system 
used in this study is characterised in C 1 23 . The purpose of this study is the 
spectral analysis of the SVF signal where the multiple reflections may be 
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overlapped Homomorphic signal processing techniques ere applied to estrect 
information about the magnitude spectrum of an individual reflection. The 
measurements on the magnitude spectra of individual reflections may be easily 
qu».ntifiably examined for NDE applications. The results of this study should 
provide more efficient and revealing CANDE schemes utilising the SVF and 
advanced digital signal processing techniques. 

Although the SVF has been generally used to characterise microstructural 
defect states of materials, it may also be used to detect cracks and 
delaminations . The output in the SVF configuration may be significantly 
affected by cracks oriented perpendicular to the faces of the test structure, 
which might be otherwise difficult to detect by many NDE methods. 
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2 SWF OUTPUT SIGNAL 

Steady- state harmonic stress wave transmission characteristics of an 
isotropic elastic plate with the transmitting and receiving transducers in the 
SWF configuration are analysed in CIO] an C Cl 11 . Circular transmitting and 
receiving transducers are considered It is assumed that the transmitting 
transducer transforms an electrical voltage into a uniform normal stress at 
the top face of the plate. Ignoring the bottom boundary, the radiation into a 
half space is considered first This introduces longitudinal (P) and shear (S) 
waves into the plate. Then the reflections which introduce reflected P and S 
waves for each incident P or S wave are considered successively at the bottom 
and top faces of the plate. For example, the multiply reflected wave specified 
by the chain PPSP is constructed from the half space solution as a P wave 
which reflects as a P wave at the bottom face, then as an S wave at the top 
face, and again as a P wave at the bottom face. The total number of P's and 
3's in the chain is denoted by p and s, respectively It is assumed that the 
receiving transducer produces an electrical voltage proportional to the 
average spatially integrated normal stress over the face of the transducer due 
to an incident wave The exact solution for the frequency response of the 
plate is formulated in integral form as the superposition of the frequency 
responses due to all the possible multiply reflected waves which reflect at 
the bottom face a final time and are then incident upon the receiving 
transducer It is found that the phases of the frequency responses due to the 
waves with common p and s have the same linear dependence on frequency Thus 
the waves with common p and s; simply the waves with p,s; have the same time 
delays in the transient analysis For example the waves with pv3,sal are the 
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the wtvn apacificd by the chains PPPS, PPSP, PSPP, and SPPP. 

The asymptotic behavior of the frequency rasponse due to a multiply 
reflected wave is also analysed. It is found that a receiving point observes 
an incident wave asymptotically as a plana wave which propagates and reflects 
in the direction of the multiply reflected ray constructed geometrically using 
Snell's law. The far field condition for the asymptotic solution to be valid 
is discussed. It suggests that although the thickness of the plate may be 
small, waves which have a sufficiently large total ray path length will 
satisfy the far field condition. Formulae for the time delays for waves 
travelling the ray path length are also given. 

The analysis given in CIO] and Cl 13 than suggests that the SVF signal 
can be modelled as the superposition of multiple reflections given by 

p + ■ -o« , s ■*« 

y (n) - y y <n-n ) < 1 ) 

P < • p » s 

pts*2,(c0 

where p*s *2 , 4 , 6 , . . , and s-0,1,2,... . Further, y(n) is the discrete- time 

5<J? signal, y (n-n > is tha sample sequence of the reflection iienaratad by 

p,s p,s r 1 

the waves with p,s; or stated otherwise, simply the reflection with p,s. n is 
an integer which defines the discrete- time location t by 

t = n T < 2 ) 

where T is the sampling period The reciprocal of T is called the sampling 

frequency n is the delay for tha reflection with p,s. Only even values of 
p.s 

p+s are considered in eqn (1) because the corresponding waves satisfy the 


r 
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condition that they reflect at tha bottom face last tiaa. 

It ia assumad that tha test spaciman and all tha componants of tha 
•xperimental ayatam ara linaar and tima invariant. Than, tha following can bo 
written C33 

y (n-n ) - x(n) * h (n) * h (n-n ) (3) 

p , s p.s a p.a p > a 

whara * donotaa tha linear convolution operation C33 . x(n) ia tha input, aignal 

r which generates tha SWF aignal h (n) ia tha sample aaquanca of tha impulaa 

* 

response of tha axporimental system without a taat apaciman . h (n-n ) ia 

P ' a p.a 

tha sample aaquanca of tha impulaa response due to tha waves with p,s. As 
discussed above, tha continuous frequency response due to a multiply reflected 
wave in tha plate ia studied in Cl 03 and Clll. Assuming that tha sampling 

frequency is greater than tha Nyquiat frequency, the FFT of h (n-n > is 

p.s p,s 

the sample sequence of tha continuous frequency response due to tha waves with 
p.s scaled by 1/T 13,121 
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3 HOMOMORP HIC PROCESSING OF SWF SIGNAL 

If there exists * signet thet is the sum of v.wo component signets whose 
Fourier trensforms occupy different frequency bends, then it is possible to 
seperete the two components by e lineer filter. On the other hend if there 
exists e signet thet is the multiplicetion or convolution of two component 
signets, then it mey be possible to seperete the two components by e 
homomorphic filter C33 

Below, the properties of e short- time portion of e SWF signet, which is 
celled the short- time SWF signel, ere discussed first. Then, the properties 
of homomorphic systems ere summenxed end homomorphic signel processing is 
epplied to the short- time SWF signel to obtein estimetes of the megnitude 
spectre of mdividuel reflections despite the feet thet the successive 
reflectionc mey be overlepped 

2_1 Short-Time SWF Signel 

For the purpose of the enelysis in this study, the complete SWF signel 

is multiplied by e window w(n) which gives the short- time SWF signel y (n> es 

w 

y(n)=w(n)y<n) ( <1 > 

w 

where w(n> is celled the dete window 

If the significent reflections in the SWF signel ere not overlepped, e 
rectenguler dete window is used end its time mtervel is chosen such thet the 


short- time SWF signel conteins only en tndividuel reflection end excludes ell 
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the ether reflections completely. Than, by taking the FFT of tihort- time SWF 
eignele the magnitude spectra of individual reflections can be found directly 
as in conventional ultrasonic spectral techniques CS3 . This is also called the 
short- time Fourier analysis Cl 33 

If the reflections in the SWF signal are overlapped, it is no longer 
possible to obtain spectral information on individual reflections directly by 
the short- time Fourier analysis However, it may be possible to obtain 
estimates of the magnitude spectra of individual reflections by the short- 
time homomorphic analysis. For this purpose the short- time SWF signal must bs 
defined appropriately as described below 

First, it is assumed that the reflections whose delays are inside the 
data window are completely covered by V w <n>, otherwise they are completely 
outside the window It is impossible to satisfy this condition esactly when 
the reflections are overlapped. However it may be satisfied approsimately by 
using a data window which is long enough and tapered smoothly to tero at the 
beginning and the end. Standard Hammmq, n -inning, or Blackman windows C33 may 
be appropriate for this purpose Then using eqns (1) and (4) 

J 

y w < n ) ~ w< n > ^ ] y.tn-n^) (5) 

j ■ 1 

where the index j replaces the complex index (p.s) in eqn (1) such that 

n =n and n,<n„<n„< <n , In this notation, the reflections in the dita 
j p , s 12 3 J 

window are ordered according to their deUys. J is the total number of 
reflections in the window y.(n-n) is eaual to y (n-n ) accordingly 


Second, it is assumed that the individual reflections in the window have 


approximately tho laat waveform but they differ only by a scale. Suppose the 
scale of the reflection y^tn-n^) Is maximum compared with the scales of the 
other reflections in the window and that it is unity Then the other 
reflections can be written as 

y^fn-n^l * A^ y^(n-n^) <*> 

where A^ is the scale of the reflection having the delay n^ and it is smaller 
than or equal to unity, given that A is equal to unity. 

Bk 

The asymptotic ray analysis in CIO! and Clll suggests that each 
reflection is generated by the stress waves in the plate which propagate and 
reflect in a particular multiply reflected ray direction, which is 
constructed geometrically via Snell's law. The direction is defined by the ray 
angle; that is, the angle between the ray and the transmitting transduisr 
axis . The frequency components of the input signal are transmitted in the ray 
directions according to the directivity functions which are given and 
discussed in detail in CIO] and Clll 

For the P waves, the increasingly higher frequency components are 
transmitteo at decreasing ray angles. The individual reflections with larger 
numbers of p, namely with larger delays, have smaller ray angles. Thus, for P 
waves, the reflections with larger delays tend to contain higher frequency 

l 

components So, v he spectral contents of two reflections may differ 
significantly if the difference between their ray angles is large; namely, the 
difference between their delays is large Thus, these two hypothetical 
reflections also differ in the time domain significantly The condition in 
eqn (6) i r then satisfied better when the length of the data window is net 
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vary long Alio, using a window which ii taparad smoothly to iaro at tha 
baginning and tha and dacraaiaa tha contribution oZ tha reflections at tha 
bag nning and and of tha short- tima SWF signal 

So , tha conditions givan in aons . (S> and (41 both ara satisf lad battar 
whan tha data window is taparad smoothly to t*.;o at tha baglntting and tha and. 
Howavar, aqn. (3) requires a lonqar window, wharaas aqn. <6) raquiras a 
shortar on# Thus, thar* is a compromise in tha langth of tha data window. 
Substituting aqn (6) into eqn. (3) 


v ( n) = win) 

w 


7) A i y m <n ‘ n i 


1 - 1 


Eqn (7) can also ba written as 


< 7 ) 


whara 


and 


y in) * win) Cv in) * pin)] 

w m 


i 8 > 


J 

pin) - ^ 

1 - 1 


A ^ <<n-n . ) 


t 9 ) 


| 1 

6 ( n-n ) ■» J 

( o 


n«n^ 
ny n 


1 


i i o > 


whara din) is called tha impulse sequence and pin) is called an impulsa train. 
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An impulsa train is dafinad aa a saquanca in which tha nonsaro aamp’aa ara 
apacad at intarvala that ara graatar than ona C143 . 

Assuming that w(n> varias slowly with raspact to y (n> , aqn . (8) can ba 

fll 

writtan as C3,131 

Y(n)*p(n)*y(n> (11) 

w w m 


whara 


p(n)«w(n)p(n> (12) 

w 

Using aqn> (9), (10) and (12), 

J 

p (n) ■ B. d(n-n) (13) 

w tmU j j 

j = l 


whara 

B. . Aj w(n.) (14) 

It is assumad (hat thu data window is locatad such that 

B . w(n ) ' 1 ( 15 ) 

m ra 

Eqn (13) is satisfiad whan tha middla of tha window is locatad naar 
tha dalay n tha casa of Hammina, Hanning, or Blackman windows C31 . 

IT 

By taking tha FFT of tha short- tima SVF signal y <n> , aqn. (11) bacomaa 

w 
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Y (k) « P (k) Y (k) (14) 

w w m 

where Y (k) , P (k) , and Y <k> are the FFT's of y (n), p (n) , and y (n>, 
w w n w w m 

respectively, k is an integer and defines the discrete- frequency location f 
by 


k 1 

f = (17) 

N f T 

where N ( is the number of points for the FFT operation, which is some power of 

two It is assumed that is adequately large so that the multiplication in 

eqn (16) corresponds to the linear convolution in eqn. (11), the conditions 

for which are discussed in C3] and [12] . 

By taking the magnitude of Y (k) , eqn (16) becomes 

w 

IY(k)l~IP(k)IIY(k)l (18) 

w w m 

where I I denotes the "maqnitude of" 

So, the magnitude spectrum of the short- time SWF signal I Y^(k) I is 
approximately equal to the multiplication of two component magnitude spectra, 
namely, the magnitude spectrum of an impulse train I P^(k) I and the magnitude 
spectrum of the particular reflection with the delay n I Y (k) I . Below, 

TO TO 

I Y^(k) I is obtained from I Y^tk) I by homomorphic signal processing. 
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3.2 Homomorphic Signal Processing 

A detailed discussion of homomorphic signal processing is given in C52 
and Cl 31. Only a summary of the relevant points is pr vided here. 

Suppose H denotes the transformation of a system such that 

x < n ) « HCs(n) 1 ( 1?) 


where s(n> and s(n) are the input and output of the system, respectively. 
Suppose that H satisfies the generalised principle of superposition under some 
rules of combinations (for example, addition, multiplication, convolution) in 
the inputs and outputs as 

HCSj(n) □ Sj(n)3 * HCs^(n)3 o HC * 2 < n) ] 

*«.<n)o* 2 <n) (20) 

where s^(n) and s^in) are two arbitrary signals. The symbols o and o designate 
the rules for combining the inputs and outputs, respectively, an^ they are 
both commutative and associative C33 . The symbols □ and o are also called the 
input and output operations, respectively. The class of systems satisfying the 
generalized principle of superposition is called homomorphic systems CIS! . 
Linear systems are a special case for which the input and output operations 
are both addition. 

Any homomorphic system can bs decomposed into three cascaded systems via 

the ordered transformations D , L, and D -1 C3,153. For the system H, the 

□ o 

system 0 Q is a homomo-phic system and has the property that 
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D [j (n) □ s,(n)] = D_C«. (n) 1 * D ts„(n>] 

□ 1 2 □ 1 □ 2 

* *j(n) + *2< n > (21) 

As observed in eqn. (21) the system D q transforms the combinations according 
to the rule □ to the additive combinations 

The system L is a linear system and has the property that 

LCs 1 (n) ♦ s j ( n ) 1 = Lts^n)] ♦ LCs^n)] 

= q^n) ♦ q j ( n ) (22) 

The system D -l is a homomorphic system and it has the property that 
o 

D ‘‘[q, (n) ♦ 3, <n>3 = D '‘tq.fn)] o D ' 1 1 q„ (n>] 

0 1 Z 0 1 o Z 

=x j (n)o* 2 (n) (23) 

As observed in eqn. (23) D 1 transforms the additive combinations to the 

o 

combinations acco ding to the rule o 

If there is a signal that is the combination of two component signals 
under the rule a . then it may be possible to separate the two components by a 
homomorphic filter whose input and output operations are both □ This 
essentially reduces to the problem of designing the linear system L C33 

Homomorphic signal processing has been successfully applied to the 
analysis of audio signals 116], images Cl 63 , speech signals 116,17,18,19], 
seismic signals 120], and echo removal 114] In this study it is applied to 
che analysis of the SWF signal 
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The schematic of tho homomorphic signal processing applied to the SVF 
output signal to estimate the- magnitude spectrum of a particular reflection is 
shown in Fig. 1. The operations applied to the SWF signal in Sec. 3.1 are also 
included in Fig . 1 . The magnitude spectrum of the short- time SVF signal 
I Y^(k) I is the input to the homomorphic system. Since I Y^flc) I is the 
multiplication of two components tP^ik)! and IY (k) I as given in aqn. (18), 
the input and output operations of the homomorphic system are both 
multiplication. Then, the corresponding systems D and D are designated as D. 
and D.~*, respectively, where the subscript denotes the multiplication 
operation. The suitable systems D., L, and D. -1 satisfying the conditions in 
eqns (19), (20), and (21), respectively, are indicated in Fig. 1. 

If the logarithm of both sides of eqn. (18) is taken 

loglY <k)l » loglP (k)l ♦ loglY (k)l (24) 

W W B 

The logarithmic magnitude spectrum (log spectrum) of an impulse train is 
rapidly varying. If the impulses are equally spaced then it is periodic C213 . 

It is assumed that the log spectrum of the reflection with the delay n , 
log I Y (k) I , is slowly varying compared with the log spectrum of the impulse 
train, loglP^(k)l . Then it is possible to separata these two superposed 
components by linear filtering C33 . Since the filtering is done on a spectrum, 
it is called frequency- invariant linear filtering and it is esactly parallel 
to the usual low- pass or high- pass time- invariant linear filtering which 
extracts either the respective slowly varying cr the rapidly varying component 
cf a signal ?n the time domain C3J 


The suitable frequency- invariant linear filtering is performed first by 
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taking the invars* fast Fourier transform (IFFT) of loglY^(k)l . Ths IFFT of 
the log spectrum is called the cepstrum C31, a term first introduced in C22] . 

So, aqn . (24) becomes 

y < n) ~ p (n) * y (n) (25) 

w w m 

where y (n), p (n), and y (n) are the IFFT's of log I Y (k) I , log IP (k) I , and 
w w m w w 

log I Y (k) I , respectively; or the cepstrums of y (n) , p (n) , and y (n) , 
m w w m 

respectively 

Since IoglY^(k)l is slowly varying and loglP^(k)l is rapidly varying, 

their IFFT's y^tn) and p^ln) occupy the low- time and high- time portions of 

the cepstrum ? w <n>, respectively. Thus, the cepstrum of y^tn) may be recovered 

from the cepstrum of ? (n) as 

w 

y (n) ~ w (n) y (n) (26) 

m c w 

where w^(n) is called the cepstrum window and covers the low- time portion of 

V (n) 
w 

The properties of the cepstrum are discussed in C33 . In particular it is 
symmetric about zero The cepstrum of a signal similar to y (n> , which has a 

Si 

slowly varying log spectrum, is concentrated near sero time and decays rapidly 

away from zero. Specifically, if y (n) is N samples long, then its cepstrum 

m y 

v^ln) occupies mainly the interval approximately from -N^/3 to hW3 C201 . 

Thus, the cepstrum window is chosen such that it covers this interval. 

The cepstrum of an impulse train is discussed in C33 , Cl 33 , and C203 If 

the impulse train is equally spaced with a spacing interval N , then its 

P 
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capstrum shows pastes at tha locations of jN whara j-± (1,2,3, . . .), which 

P 

raflacts tha pariodicity of tha log spactrum of tha impulsa train. So, y (n> 

and p (n) ara wall saparatad in y (n) whan N is graatar than N /3, which 
w *w p * y 

limits tha validity of aqn. (26). 

If tha impulsas ara locatad arbitrarily and thay hava amplitudas of 
arbitrary strangth and random polarity, tha impulsa train may ba thought of as 
a noisaliks random sigr.^l. Than its log spactrum is rapidly varying but not 
pariodic . Thus its capstrum doas not show any paaks. Howavar, aqn. (26) may ba 
still valid C3.13V 

Ones y (n) Is racovarad by aqn. (26), loglY (k) I is found by ta v ing tha 

A A 

FFT of y (n) , and than I Y (k) I is found by ar.oonantiation of log I Y (k) I . 

A A IQ 

log I Y (k) I , which is tha log spactrum of tha raflaction y (n) , is also callad 
in a 

tha capstrally smoothad log spactrum of tha short- tima SWF signal y (n) . 

w 

Tha tima history of tha raflaction y (n) cannot ba found by tha procass 

IQ 

shown in Fig . 1 bacausa tha logarithm oparation is parformad on I Y <k) I and 

w 

the phasa information is lost . If tha complas logarithm oparation is parformad 

on Y^(k) than it may ba possibla to find y^ln) by adding anothar IFFT 

oparation to tha and of tha procass. Tha complaa logarithm raquiras unwrapping 

of tha phasa of Y (k) to obtain a continuous phasa curva C33 . Tha 

w 

corresponding term for tha capstrum is callad tha complas capstrum, which is 
the IFFT of tha complas logarithm of tha FFT of a tima domain signal. For tha 
analysis in this study, tha capstrum is usad rathar than tha complas capstrum. 


A 
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A RESULTS AND DISCUSSIONS 

The experiment*! system used in this study is characterised in C123 by 
coupling the transmitting and receiving transducers face- to- face, without a 
test specimen in between. The impulse response of the experimental system and 
its loq spectrum are shown in Figs. 2 (a) and 2 <b> , respectively. The signal 
in Fig. 2 (a) and all the time domain signals in the subsequent figures are 

discretes their sampling period is 0.03 ksec (sampling frequency 20 MHt) . The 
log spectrum in Fig. 2(b) and all the frequency domain plots in the subsequent 
figures are also discrete. They are generated by a 512-point FFT. The points 
in all the time and frequency domain plots are connected by linear 
interpolation . 

The time history of the input signal to generate the following SWF 
signals, and its log spectrum are shown in Figs. 3 (a) and 3 (b), 
respectively Two tests are considered as typical examples for the SWF 
signal In Test I, the thickness of the plate h is 2.4 cm and the distance 
between the transmitting and receiving transducer axes < is 6 cm. In Test II, 
h is 1.25 cm and f is i. 25 cm. The plates are made of aluminum (6061-T6), 
which has P wave and S wave velocites of 0 . 432 and 0.J13 cm/ p sec, 
respectively The transducers (Panametrics model V105) have a radius of 1.125 
cm The SVF signals for Tests I and II are shown in Figs. 4 (a) and 4 (b) , 
respectively The time delays for some possible reflections are calculated by 
the formulae given in Cl 1 3 considering the multiply reflected rays connecting 
the centers of the transducers, and their locations are indicated above the 
time signals with the values of p and s 

The ray angles, total ray path lengths, and time delays for some 
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reflections with only P wavti (««0) in Toots I and II art given in Table 1 . 
These are calculated using the formulae given in CIO! and Clll. The distances 
between the transducer aies are chosen such that the ray angle for a 
reflection with pap^.saO in Test I is the same as the ray angle for the 
reflection with pa2pj,s-0 in Test II, where p^ is an arbitrary positive even 
integer As observed from Table 1 for some selected values of p , their total 
ray path lengths and time delays are also approximately the same. Then, it is 
expected from the analysis in CIO! and Cl 13 that these two corresponding 
reflections have the same waveform because all the frequency components of the 
input signal are transmitted to the receiving transducer in the same direction 
as discussed in Sec. 3 

The reflection with p«Pj , s*0 in Test I and the reflection with p^p^ , 
s*0 in Test II are generated from the stress waves which reflect successively 
at the bottom and top faces of the plate <p^-l) and <2p^ — 1> times,, 
respectively. Thus, it is also expected that the reflection with pa2p^ , saO in 
Test II is approximately (Qpp)^ times smaller than the reflection with p»p^ , 
s = 0 in Test I, where Qp p is the reflection coefficient for a P wave to a P 
wave reflection. 

There are several observations which can be made from Figs. 4 (a) and 4 
(b) In Test I, although there are many possible P and S wave multiple 
reflections, only the multiple reflections of the P waves (s»0) and of the P 
waves containing a single S wave (sal) are above the noise level. Most of 
these reflections are not overlapped significantly . This nonoverlapping case 
exists because the thickness of the plate is sufficiently large that the 
difference between the time delays for the successive significant reflections 
is adequately large Thus the short- time Fourier analysis is suitable for 


1 


20 


obtaining spectral information on individual redaction* in the SWF signal 
shown in Fig . 4 (a) . 

In the SWF signal for Test II shown in Fig . 4 (b) , successive 

reflections are overlapped and thus cannot be identified separately. This is 
because the thickness of the plate is so small that the difference between the 
time delays of the successive significant reflections is smaller than the 
length of a single reflection. Here, the short- time Fourier analysis does not 
give a direct access to spectral information. However, as discussed in See. 4, 
the short- time homomorphic analysis gives an estimate of the magnitude 
spectrum of a single reflection although it is overlapped with other adjacent 
(and perhaps beyond) reflections. 

The short- time Fourier analysis of the SWF signal for Test I is 
obtained and spectral information on the reflection with p-12, s-0 is shown 
in Fig. 5 The time history of the short- time SWF signal, which is the 
reflection with p>12, s = 0, and its log spectrum are shown in Figs. S (a) and 5 
(b) , respectively The characteristics of the data window are indicated in the 
title for Fig. 5 (a). The origin of the abscissa (time) in Fig. 9 (a) is 

shifted to the beginning of the data window 

The short- time homomorphic analysis of the SWF signal for Test II is 
shown in Fig. 6. Here the goal is to obtain spectral information on the 
reflection with p*24, s*0 The process in Fig. 1, which is discussed in Sec. 
4, is used The short- time SVF signal and its log spectrum, cepstrum and 
cepstrally smoothed log spectrum are shown in Figs 6 (a), (b>, (c), and (d) , 
respectively Considering the discussions in Sac. 3, a Hamming window is used 
as the data window, and its beginning and end are chosen such that it covers 


the reflection with p-24, s = 0 in its middle part. Since there is a compromise 
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in selecting the length of the window to satisfy the conditions for eqns. (5) 
end (6) simultaneously , the length of the date window is chosen by trial. 

The log spectrum in Fig. 6 (a) contains rapid variations superposed onto 
what appears to be a slowly varying component. As discussed in Sec. 4, the 
rapid variations come from the impulse train which contains impulses located 
at the delays of the reflections in the short- time SWF signal. 

Notice that the cepstrum in Fig. 6 (c) is substantially confined near 
zero It also contains small peaks at ± 3.8 Msec away from sero. As discussed 
in Sec. 4, the small peaks away from zero are due to the impulse train and 3.8 
Msec is approximately the same as the difference between the time delays of 
the reflections with only P waves. 

The cepstrum in Fig. 6 (c) is multiplied by a Hamming window having a 
duration from -1.25 Msec to 1 25 Msec, which is called the cepstrum window, 
and then its FFT is taken to obtain the smoothed log spectrum in Fig. 6 (d) . 
The Hamming window is chosen by trial. Considering the discussions in Sec. 4, 
its half length, 1 25 Msec, is chosen approximately one third of the length of 
tne reflection with p«12, s*0 in Test I. The cepstrally smoothed log spectrum 
in Fig 6 (d) is an estimate of the log spectrum of the reflection with p»24, 
s = 0 in Test II. The log spectrum in Fig. 6(b) and the cepstrally smoothed log 
spectrum in Fig. 6 <d) are superposed in Fig. 6(e). 

The short- time Fourier analysis of the SWF signal for Test I is 
repeated with different data windows to obtain log spectra of the reflections 
with p*8, s = 0; and with p = 10, s = 0 Also, the short- time homomorphic analysis 
of the S«/F signal in Test II is repeated with different data windows to obtain 
estimates to the log spectra of the reflections with p-16, s*0; and with p*20, 
s-0 
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The characteristic* of the log spectra of all tha raflsctions considered 
above are listed in Table 2. The characteristics of the data windows are also 
indicated. The parameters used in Table 2 such as the peak log magnitude, peak 
frequency, average resonant frequency, lower and upper frequencies of -3 dB 
and -4 dB frequency bands are defined in Fig. 7 on a hyp<>thatical log 
spectrum 

As discussed at the beginning of this section, since the reflee’.ion with 
P^P^, s-0 in Test 1 and the reflection with p«2p^ , s-0 in Test II have the 

same ray angles, their log spectra are espected to have the same shapes but 
different scales The characteristics of these two corresponding reflections 
are compared in Table 2 with the relative errors given in parentheses. The 
definition of the relative error is given in the Table 2 footnotes. 

Notice that the -3 dB and the -4 dB frequency bands of the reflections 
in Tests I and II tend to move to higher frequencies for increasing values of 
p As discussed in Sec. 3, this is because the reflections with larger values 
of p have smaller ray angles and tne higher frequency components of the input 
signal are transmitted to the receiving transducer at smaller ray angles 
according to the directivity functions discussed in C13] and Cl 13 . 

The discussions at the beginning of this section also suggest that the 
reflection with p-2p^ , s*0 in Test II is espected to be (Qpp)** times smaller 
than the reflection with p-p^, s»0 in Test I. The values of 20 log(Qpp> P are 
read from Table 2 as the differences between the peak log magnitudes of the 
corresponding reflections in Tests II and I Then, an esperimental reflection 
coefficient for P waves to P waves is calculated for each value of p^ 
considered . 


I 

1 
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The esperimental reflection coefficients are compared with the 
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theoretical reflection coefficient* in Tab!* 3. The thaoratieal reflection 
coefficient! are calculated fro* the formulae given in CIO] and C231 for the 
reflection of plane P waver at a i*.re«- free plana boundary. The incidence 
angle used in the formula ie the same as the ray angle given in Table 1 for 
each value of p^ considered. Notice that the esperimantal reflection 
coefficients are smaller than the corresponding theoretical reflection 
coefficients This is consistent with the espectatlon that there is likely to 
be an energy loss (apart from mode conversions) in reflections. 

Finally, two hoik tests are considered to demonstrate how the short- 
time homomorphic analysis given in this study may be utilised in the NDE of 
plates containing cracks perpendicular to its surfaces. For Test III ,the 
plate is (he same as for Tost II (h«l 25 cm, aluminum) and < is A cm. For Test 
IV, the same slate, but now containing a crack, is used, and < is again 6 cm. 
The crack is open to the bottom face of the plate like a slit all the way 
along the depth (into the plane of the sketch) of the plate. It has a depth d 
of 0 62 cm and a width of 0 05 cm The crack faces are perpendicular to the 
crack and is located midway between the transducers. 

The SWF signals for Tests III and IV are shown in Figs. 8 (a) and I (b), 
respectively. The time history of the short- time SWF signals, their log 
spectra and their cepstrally smoothed log spectra aie shown in Figs. 9 (a) and 
9 (b) for Tests III and IV, respectively Considering the location of the data 
window, the cepstrally smoothed log spectra represent the estimates of the log 
spectra of the reflections with p-18, s>0 in both Tests III and IV. These two 
log spectra are compared in Table 4 The relative differences between the 
defined characteristics of the two log spectra are given in parentheses. The 
relative difference is defined under the Table 4 footnotes No significant 
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reletie* differences in these defined ehereeterieties ere observed ee 
indicetors of the creek. However, the peek log aegritudes in Tests III end 
IV ere -1.03 end -0.33. respectively. So, the peek log aegnitude drops 7.32 
dB (the peek aegnitude drops to 42%) due to the existence of the treck. 
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The SWF output signal was analyzed in the frequency domain. The SVF 
signal, which is the superposition of successive multiple reflections, was 
multiplied by a dat.i window and its short- time portions were analyzed i« 
obtain spectral information contained in individual reflections. 

If the significant reflections were not overl pped, the short- time 
Fourier analysis was used, that is, the short- time SVF signal contained only 
a single reflection and its FFT was computed . 

If the reflections were overlapped, it was assumed that the reflections 
in the short- time SVF signal have approximately the same waveform but 
different scales This assumption suggested a data window which was tapered 
smoothly to zero at beginning and the end of the window and long enough to 
cover the included reflections completely, but not so long that the waveforms 
of the included reflections changed significantly. Then, the short- time SVF 
siqnai was modelled as the convolution of an impulse train, where the impulses 
were located at the time delays of the reflections, and the reflection with 
the maximum scale 

A summary of the relevant points of homomorphic signal processing was 
given Homomorphic signal processing was applied to short- time SVF signals to 
obtain estimates of the log spectra of individual reflections although the 
individual reflections were overlapped 

Two typical SVF signals generated in aluminum plates were analysed. In 
Test I, the SVF signal contained nonoverlapping reflections In Test II, the 
thickness of the plate was approximately half of that in Test I, and the 


reflections were overlapped The reflection which was generated from the 
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stress wivti which travelled p times as a P wave and s times as an S wave 

through the plate while reflecting back and forth between thu. bottom and top 
faces of the plate was designated as the reflection with p,s. The multiply 
reflected rays which were constructed geometrically between tV centers of the 
transmitting and receiving transducers using Snell'? law were considered The 
distances between the transducer ases in Tests I and II were chosen such that 
the reflections with p-p^ *«o in Test I and the reflections with p-2p t . e-0 

in Test II had the same ray angles and approsimately the same total ray path 
lengths The theoretical analysis suggested that the magnitude spectra of the 
two corresponding reflections with the same ray angles in Tests I and II were 
expected to have the same shape but different scales because all the frequency 
components of the input signal are transmitted to the reflections from the 
same directions. 

The log spectra of the reflections with p-p^ , s«0 in Test I were 
obtained by the short- time Fourier analysis for the values of p^ of 8, 10, 
and 12 The estimates of the log spectra of the corresponding reflections with 
P*2Pj , s = 0 in Test II were obtained i >7 the short- time homomorphic analysis. 
The characteristics of the log spectra the two corresponding reflections in 
Tests I and II were compared in terms of their peak log magnitude^, peak 
frequencies, average resonant frequencies, and lower and upper frequencies of 
-3 dB and -6 dB frequency bands It was observed that in general they were 
approximately the same except for the peak log magnitudes, as the theory of 
wa/e propogation in isotropic elastic plates suggested. 

From the differences between the peak log magnitudes and between the 
values of p in Tests I and II, experimental reflection coefficients were 
calculated for different ray angles These showed good agreement with the 
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theoretical raflaction coefficients for plana wavas. Alao, it was observed 
that fraquancy bands of tha raflactions with largar p, namely with smaller ray 
anglas, moved to highar frequencies, which was consistant with tha bahavior of 
tha diractivity function of tha P wava. 

Finally, tha potantial usa of tha results of this study for tha NDE of 
plates containing cracks perpendicular to their faces was demonstrated in a 
specific test. For a specific set of parameters such as tha input signal, 
plate thickness, crack depth, location, and orientation, it was observed that 
tha decrease in the peak log magnitude of a specific raflaction, which was 
found by tha short- time homomorphic analysis, was tha indicator of tha crack. 

This study provides a quantitative analysis of the SWF signal in tha 
frequency domain . Since the information in tha frequency domain can be related 
directly to tha steady- state wava propagation analysis in test structures, it 
should provide more revealing measurements for NDE applications utilising the 
SWF The effect of different cracks on the characteristics of tha log spectra 
of different reflections should be analysed. This study may be useful in the 
analysis of the ultrasonic output signals in different test structures which 
may contain overlapped reflections. Ultrasonic attenuation measurements and 
transducer calibration techniques may benefit from this study. The short- time 
homomorphic analysis may also be applied to acoustic emission (AE) signals to 
characterise unknown sources in terms of their directivity functions Finally, 
developments in CANDE should be faciliated by the digital signal processing 
procedures which are summarised in this study 
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Soma Raflactiona With Only P waaas in Taata I and II . 
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30 44 


48 44 

















TADLE 2 ChtracltucUcs of Log Bpoctio of Soma Rofloction* With Only P Wivos 
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TABLE 3 Comparison of Experimental and Theoretical Reflection 
Cosff lelsnts . 


Valus of p 

"*T 

Angle 

(Osgrss) 

Theoretical 

Rsflsetlon 

Costflcisnt 4 

< ^PF > t 

Espsrlmsntal 

Rsflsetlon 

Cosffielsnt 

(( V. 

Rslativs 
Error b 

(%> 

Tost I 

Tsst II 

• 

14 

17 33 

0 713 

0 .144 

7.B 

10 

20 

14.04 

0 743 

0 170 

3 4 

12 

24 

It 77 

0.740 

0.703 

3 7 


a Absoluts values are oiven 
b Relative Error » l (Q pp> t “ <Qpp ) > 1 1 1 ^ppV 
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Fig S Short-time Fourier analysis of SVF output signal for Toot I 
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16 Abstract 

The stress wave factor (SWF) signal, which is the output of an ultrasonic testing system, where the 
transmitting and receiving transducers are coupled to the same face of the test structure, is analyzed 
in the frequency domain. The SWF signal generated in an isotropic elastic plate is modelled as the 
superposition of successive reflections. The reflection which is generated by the stress waves which 
travel p times as a longitudinal (P) wave and s times as a shear (5) wave through the plate while re- 
flecting back and forth between the bottom and top faces of the plate is designated as the reflection 
with p, s. Short-time portions of the SWF signal are considered for obtaining spectral information on 
individual reflections. If the significant reflections are not overlapped, the short-time Fourier 
analysis is used. A summary of the relevant points of homomorphic signal processing, which is also 
called cepstrum analysis, is given. Homomorphic signal processing is applied to short-time SWF sig- 
nals to obtain estimates of the log spectra of individual reflections for cases in which the reflec- 
tions are overlapped. Two typical SWF signals generated in aluminum plates ar analyzed. In Test I, 
the SWF signal contains nonoverlapping reflections. In Test II, the thicknes' of the plate is approx- 
imately half that in Test I, and the reflections are overlapped. Multiple rei ections of only P waves 
are considered. It is observed that the log spectra of two corresponding reflections witn the same 
ray angles and approximately the same total ray path lengths in Tests I and II have approximately the 
same characteristics except for their overall scales, as the theory of wave propagation in isotropic 
elastic plates suggests. From the differences between the overall scales and the values of p of the 
two corresponding reflections in Tests I and II, experimental reflection coefficients for P wave to P 
wave reflections are obtained for the different ray angles considered. They show good agreement with 
the theoretical reflection coefficients for plane P waves. Also, it is observed that the frequency 
response at the maximum value of the log spectrum moves toward nigher frequencies in the reflections 
with larger p, namely with smaller ray angles. This is consistent with the behavior of the directiv- 
ity function. The potential use of these results for the nondestructive evaluation of plates contain- 
ing cracks perpendicular to their faces is demonstrated in a specific test. 
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